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We investigate a supersolid state in hardcore boson models on the face-centered-cubic (FCC) 
lattice. The supersolid state is characterized by a coexistence of crystalline order and superfluidity. 
Using a quantum Monte Carlo method based on the directed-loop algorithm, we calculate static 
structure factors and superfluid density at finite temperature, from which we obtain the phase 
diagram. The supersolid phase exists at intermediate fillings between a three-quarter-filled solid 
phase and a half-filled solid phase. We also discuss the mechanism of the supersolid state on the 
FCC lattice. 
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Whether a supersolid state, where both of solidity and 
superfluidity coexists, realizes in matters or not. The 
possibility of the supersolid in ''He have been discussed 
theoretically[l|, [2, H, 0) Hi- Recently, characteristic be- 
haviors of superfluidity have been reported in torsional 
oscillator experiments on the solid ^He by E. Kim and 
M. H. W. Chan(KC)@, while a number of experiments 
failed to detect the superfluidity in the solid ''He. In their 
measurements, a sudden drop in the resonant period was 
observed around T ~ 0.2K. The drop implies emergence 
of non-classical rotational inertia in the solid ''He. They 
concluded that this was a signature of a transition into 
the supersolid phase. However, some experimental and 
theoretical studies suggested different interpretations of 
KC's observation. For instance, the signal of the super- 
fluidity in the solid ^He became weaker as annealing cy- 
cles were repeated*^ and the superflow was blocked by 
the solid ^He with no grain boundaries [SJ. Theoretically, 
the possibility of the superflow induced by vacancies in 
the commensurate solid ^He, in which the total num- 
ber of atoms equals a multiple of the number of lattice 
sites, was ruled out[l, and the superglass behavior 
appeared in a quenched system[ll|. Although several 
mechanisms of the superflow in the solid '^He have been 
proposed, the satisfactory interpretation of KC's results 
is still controversial. 

Bosonic lattice model was introduced as a reason- 
able model of liquid ''He[l2|- Recently, the possibil- 
ity of the supersolid on the lattice model in triangu- 
lar lattice [H, [3, [H, [ll] and kagome lattice cases 
was studied by quantum Monte Carlo simulations. From 
these studies it became clear that the frustrated interac- 
tions on the triangular lattice stabilize the super current 
induced by vacancies in the crystalline ordering with the 
wave vector Q=(47r/3,0) or (27r/3, 0). In contrast, the 
supersolid is not stabilized on the kagome lattice where 
the frustrated interactions exist. In the three dimen- 
sional lattice cases, the phase diagrams of the system on 
the body-centered-cubic (BCC) lattice was obtained by 
a mean-field approximation and concluded that the su- 
persolid state appears if the next- nearest- neighbor inter- 
actions are present[l,[l|. However, the reason for the sta- 
bilization of the supersolid state on the BCC lattice and 



the microscopic picture was not cleared from the mean- 
field results. (Note that the BCC lattice is bipartite and 
has no frustration if one docs not take into account the 
next-nearest-neighbor interactions.) 

Theoretical study beyond the mean-field theory for 
three dimensional systems is still missing. In this letter, 
the supersolid state in a three-dimensional bosonic-lattice 
model is studied by a quantum Monte Carlo method 
based on the directed loop algorithm [Toj . We wish 
to address a generic question what ingredient is neces- 
sary to realize the supersolidity. From the study of the 
two-dimensional case mentioned above, it is presumable 
that the geometrical frustration plays an essential role 
in the supersolidity in the bosonic lattice model. We 
therefore focus on a hardcore-bosonic model on the face- 
centered-cubic (FCC) lattice, which does not have a di- 
rect connection to the lattice structure of the real solid 
helium, because the FCC lattice is one of the simplest 
lattices with geometric frustration. 

More specifically, we consider bosonic lattice model 
with the positive hopping amplitude t > and the 
nearest-neighbor repulsion V > on the FCC lattice. 
The model Hamiltonian is defined by 
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where ji is the chemical potential, bi^{bi) is the bosonic 
creation (annihilation) operator, and fii — hi^bi. The 
summation {ij) is over the nearest neighbor pairs and the 
system size is defined by A'^ = L'^. The periodic boundary 
condition is applied. Under the hardcore condition, the 
original bosonic-lattice model is identically mapped onto 
the 5*= 1/2 XXZ model, 

(y> {ij) 
-HY^S,^, (2) 
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where J± = 2t, Jz = —V and H = p. — 6y. Note that 
and < (>)0 mean the antiferromagnetic (ferro- 
magnetic) interactions. In the spin language, the super- 
solidity is characterized by the following two properties: 
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sublattice-dependent expectation values of the longitudi- 
nal spin components (broken translational symmetry, or 
crystallization), and non- vanishing transverse spin com- 
ponents (ofF-diagonal long range order, or superfluidity). 

In the limit J±/\Jz\ — > (Ising model), the or- 
dered states were investigated and the H-T phase di- 
agram was obtainedflO, l2ll . [2^ . At the magnetiza- 
tion 'm={J2i Si^ /N)=0 and 1/4, there appear two solid 
phases conventionally referred to as AB and A3B. Repre- 
sentative spin configurations of the two phases are shown 
in Fig. [T] The phase transition from the AB phase to 
the A3B phase occurs at i?ising — 2|Jz| at absolute zero 
temperature. 




FIG. 1: The perfectly ordered spin configurations, (a) AB 
state and (b) A3B state. Note that the dashed lines connect- 
ing next-nearest neighbors are mere guide lines to the eye and 
there is no direct coupling corresponding to these lines. (The 
direct couplings exist only for nearest neighbor pairs.) The 
gray lines denote the superflow paths in the supersolid state 
(see text). 



In order to investigate the crystalline order and the 
off-diagonal long-range order for J± > 0, we calculate 
the static structure factor S{Q) (SSF) and the superfluid 
density ps [Hi , defined by 



SiQ) = 



and 
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where W = {Wx,Wy,Wz) denotes the winding number 
of the world-lines. In what follows, we express the wave 
vector Q by the conventional choice of the unit reciprocal 
vectors. 

When the system is in the AB or A3B ordered state, 
the SSF is proportional to N and strong system-size 
dependence is expected at Qsoi=(7r, tt, 0), (7r,0, tt), and 
(0, TT, tt). For other Qs, the SSF should be system- 
size independent. It was reported for the classical case 
(Jj.=0)[22, [2I] that a perfect sohd state, either AB or 
A3B, can hardly be observed in a system of computa- 
tionally accessible size due to antiphase domain bound- 
aries (APB). While the states with domain boundaries 
have negligible weight in the thermodynamic limit, they 
have non-negligible contributions for small systems be- 
cause the domain-wall free-energy is small due to the 



frustrated nature of the interactions. The APBs reduce 
'S'(Qsoi) since contributions from different phases have 
opposite signs. However, the magnitude of the reduction 
depends on the locations of the APBs and the cancella- 
tion does not in general make 5(Qsoi) completely vanish- 
ing. Therefore, the average ^(Qsoi) is still proportional 
to the system size even if the APBs are present. As for 
the effect of the APBs on the superfluid density, we have 
confirmed that it is relatively minor compared to that on 
'S'(Qsoi)- To see this, we evaluated ps in two ways (see 
Fig. 4 (b)). One is a long equilibrium simulation start- 
ing from random initial configurations, in which APBs 
are observed. The other is relatively short Monte Carlo 
simulations starting from the perfect AB or A3B con- 
figuration. In the latter, the length of the simulation is 
chosen such that APBs do not appear. In both cases, 
the superfluid density yielded the same value within the 
statistical error. 

In Fig. [21 we show the results of the field depen- 
dence of S'(Qsoi) and ps at {J±,Jz) = (0.2,— 1.0) J and 
ksT = O.IJ, where J = jj^j is our unit of energy. 
The _ff-axis can be divided into four regions according to 
the behaviors of S'(Qsoi) and ps] (I) the low-field region 
< H < TJsoiidi ^ 1.15J, (II) the lower-intermediate 
region ffsoiidi < H < H* ^ 2.2J, (III) the upper- 
intermediate region H* < H < -ffsoiid2 ^ 3.1J, and (IV) 
the high-field region iJsoiid2 < H, where iJsoUdi, H* , and 
^?soiid2 are temperature-dependent transition fields. 

In the regions I, III and IV, the crystalline order ex- 
ists. This is evident from the fact that S'(Qsoi) increases 
in proportion to the system size. In the region II, on the 
other hand, 5'(Qsoi) does not show a system-size depen- 
dence indicating no crystalline ordering in this region. 
The superfluid density ps is almost zero in the regions I 
and IV whereas in the intermediate regions II and III, it 
stays finite. Judging from these results, we conclude that 
the ground state is the solid state in the regions I and IV. 
As shown in Fig. [3l the magnetization plateaus at m—Q 
and 1/4 appear in the corresponding fields. Therefore, 
these solid phases are the AB and A3B ordered phases, 
respectively. The region II is the superfluid phase. Fi- 
nally in the region III, since the crystalline order and the 
superfluidity coexist, there appears the supersolid phase. 
Hence, we conclude that the supersolid state is stable in 
three dimensions. 

Next we study the temperature dependence. The re- 
sults at 7? = 2.7J (in the supersolid region III) are shown 
in Fig. m As we decrease the temperature with fixed 
magnetic field, 5(Qsoi) almost discontinuously increases 
at kBT=kBTsoiid ^ 0.32J and the system-size depen- 
dence appears for T < TsoUd- However, the superfluid 
density remains very small; < 5 x 10~^ near TgoUd- 
This is the transition from the normal fluid to the solid 
phase. Since S'(Qsoi) takes the same values of those in the 
A3B solid region IV and scarcely shows the temperature 
dependence in T < Tgoiid, we identify this solid phase as 
the A3B ordered phase. At even lower temperature, the 
system undergoes another transition from the solid phase 
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FIG. 2: The field dependence of S'(Qsoi) and ps at (Jj_, Jz) = 
(0.2, — 1.0)J and fcsT = O.IJ. Tlie open circles, the solid cir- 
cles, and the inverted triangles denote the results of L = 6, 
12, and 18, respectively. Note that the results of 5'(Qsoi) 
are averaged values of S{Q,TT,-n), S{Tv,0,n) and S(7r,7r,0). 
The inset shows the system size dependence of S{Qso\)/N 
at i//Jz = 3.0. 



results are shown in Fig. [5] 
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FIG. 4: The temperature dependence of S'(Qsoi) and ps for 
{J±, Jz)=(0.2, -1.0) J and H=2.7J. In (b), the cross symbols 
denote the results starting from the perfect A3B configuration 
in L=18 and the others are those starting from the random 
initial configurations. The inset in (b) is the finite-size scaling 
of the superfluid density. 
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FIG. 3: The field dependence of the magnetization at 
(J±, Jz) = (0.2, -1.0)J and kBT=0.1J. The value of m=l/4 
corresponds to a half of the saturation magnetization. 
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to the supersolid phase. This is marked by the increase in 
the superfluid density ps that starts at fcsTsupor ^ 0.22 J. 
To estimate Tguper: we analyze the finite-size-scaling be- 
havior of the superfluid density ps by the scaling form 
PsL=f{L^/''{T — TsupcT:)) using the exponents of the three 
dimensional XY model z^=0.6723[2^. As shown in the 
inset of Fig. [IJb), the data collapse is obtained with the 
critical temperature, fcBrsuper=0.221(2) J. Thus, we con- 
clude that the phase transition from the solid phase to 
the supersolid phase is of the second order and its uni- 
versality class is that of the three dimensional XY model 
as expected. In this way, we estimate the critical tem- 
peratures Tsoiid and Tgupor for various other values of H 
to obtain the phase boundary. Here, Tgoiid and Tsupcr 
denote the transition temperatures where the crystalline 
order and the superfluid order, respectively, emerge. The 



FIG. 5: The H-T phase diagram for (Jx, Jz) = (0.2, -1.0) J. 
The open circles and the solid circles indicate the first-order 
and the second-order transitions, respectively. The solid lines 
are mere guide lines to the eye. The labels, "N.F.", "S.F.", 
and "S.S." stand for the normal fiuid phase, the superfluid 
phase, and the supersolid phase, respectively. 



From the phase diagram in FigO we flnd that the 
supersolid phase exists between the AB and A3B solid 
phases. This region locates slightly above i?ising=2.0J 
at which the phase transition occurs in the classical case 
{J±—0) from the AB to A3B phase. To discuss the mech- 
anism of the supersolid state on the FCC lattice, we con- 
sider the spin configuration above i?ising- In the classical 
case, at the critical field, the spins at the centers of faces 
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of the cubic lattice (the up spins in Fig. [Ijb)) become 
dangling spins; they can be reversed without changing 
the energy. Let us regard down (up) spins at these lo- 
cations as hard-core particles (holes). Then, H — i?ising 
can be interpreted as the excitation gap for creating a 
particle and the ground state is the empty state at H 
larger than the critical value. However, once the hop- 
ping (J_l) is turned on, the excited particles may move 
along the gray lines in Figlljb) and can in general con- 
dense. If the magnetic field is far larger than Hisingj 
the classical gap is larger than the scale of the hopping 
constant J± (i.e., the band width of particle excitation), 
the gap remains open even if the quantum hopping is 
present. As we decrease the magnetic field, however, at 
some point the classical gap becomes smaller than the 
scale of the hopping constant. Accordingly the actual 
gap closes and the ground state starts exhibiting super- 
fluidity. At this point, in contrast to the spins on the 
faces, the spins at the corners can hardly be affected by 
the hopping term, because the energy cost of reversing 
one of these spins is AE^ j\Jz\ and is still too large. 
Therefore, they stay in a solid crystalline order. As we 
further decrease the magnetic field, the density of con- 
densed particles at the dangling spin locations gradually 
increases. This generates positive molecular fields at the 
corners, destabilizing the crystalline order. This desta- 
bilizing effect finally melts the crystal at H—H*. This 
latter transition point must be larger than -ffising because 
the transition must take place before the classical exci- 
tation gap closes and therefore the density of the excited 
particles diverges. This is the microscopic scenario of the 
two transitions, the solid to the supersolid transition and 
the supersolid to the super fluid transition. Indeed, we 
successfully confirmed the existence of the AaB-type su- 
persolid states in the corresponding parameter region in 
the present simulation. 

This scenario predicts a supersolid phase of another 
type in the region H < ifising, which we could call AB- 
type. The mechanism of the AB-type supersolid is again 
understood by the dangling spins. This time, the dan- 
gling spins appear at the sites occupied by down spins 
in FiglTJa), and we should regard the up spins on these 
sites as excited particles, which condense in the AB-type 
supersolid phase that locates below i?ising- The excited 
particles hop along the gray lines in Fig[lja), while rigid 
spins (those pointed up in FiglIJa)) stay in the crys- 
talline order. The AB-type supersolid has the charac- 
teristic two-dimensional paths of the superfluid due to 
the alternatively stacks of the superfluid and solid layers, 
while the AsB-type supersolid has the three dimensional 



superfluid connections. Perturbatively, the effective in- 
teractions between these superfluid layers may arise in 
the second order of Jj_. In Fig. 6, we show some re- 
sults of the superfluid density and the structure factor at 
(Jj_, J^)=(0.15,-1.0)J and H^IMJ < i/hing- In this 
case, while an anomaly, which is cleared by the calcu- 
lations in the larger system size, appears in S'(QsOl) at 
ksT 1.5 due to the APB's, the AB-type supersolid 
realizes in the region ksT < fcsTsupor ~ 0.104(2). 
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FIG. 6; (i)The temperature dependence of ^(Qsoi) and ps in 
(Jx, Jb) = (0.15, -1.0) J and H^1.65J. (ii)The cross-section of 
a spin configuration in the AB-type supersolid at ksT — 0.1. 
The gray circles denote up spins and the unoccupied sites 
correspond to down spins. 



To summarize, we have calculated the SSF and ps for 
5"= 1/2 XXZ model on the FCC lattice and obtained a 
phase diagram at fixed J±/Jz- We have also discussed 
the microscopic mechanism of the supersolidity in the 
present model and pointed out that the connections of 
dangling spins resulting from the geometrical frustration 
play a key role in the formation of the supersolid state. 
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